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We study quantum interference of electronic waves that are scattered by multiple non-magnetic 
impurities in a d^2_j,2-wave superconductor. We show that the number of resonance states in the 
density-of-states (DOS), as well as their frequency and spatial dependence change significantly as 
the distance between the impurities or their orientation relative to the crystal lattice is varied. Since 
the latter effect arises from the momentum dependence of the superconducting gap, we argue that 
quantum interference is a novel tool to identify the symmetry of unconventional superconductors. 
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Over the last few years, the study of impurities 
in unconventional suaerepuductors has attija|eted con- 
siderable theoreticalaHu'Ll and experimentalElEfLl atten- 
tion. In particular, a series of groundbreaking scanning 
tunneling microscopy (STM) experiments has provided 
detailed information on|-|the density-oLstates (DOS) 
near single non-magneticD and magneticu impurities in 
Bi2Sr2CaCu208+5, a high-temperature superconductor 
(HTSC). Of particular interest is the experimentally ob- 
served emergence of resonance states near the impuri- 
ties. Several theoretical scenarios have been proposed 
that ascribe the origin of these resonances to electronic 
scattering off classical impuritiea-|rwith magnetic and 
non-magnetic scattering potentialsaHu or to the Kondo- 
screening of a local spinj-polarization that is induced by 
non-magnetic impuritiesQ. At the same time, a different 
line of beautiful experiments studied quantum interfer- 
ence of electronic waves that are scatteredJpy multiple im- 
purities. In particular, Manoharan et al.u demonstrated 
that quantum interference in a corral of magnetic impu- 
rities arranged on the surface of a metallic host leads to 
the focusing of electronic waves into a quantum image; a 
result that has recently been addressed in several theoret- 
ical studiesu. Similar quantum interference experiments 
using unconventional superconductors can be expected 
in the future. First evidence for quantum interference in 
the HTSC was recently reported by Derro et alB in the 
one-dimensional chains of YBa2Cu306+x. 

Motivated by the experimental progress in this field, 
we present in this article a theoretical model which com- 
bines the study of impurities in unconventional super- 
conductors with that of quantum interference effect 
particular, generalizing the formalism presented iro'ElE. 
we investigate the electronic structure in the vicinity of 
two non-magnetic impurities in a ci2,2_j,2-wave supercon- 
ductor. We show that quantum interference due to the 
presence of a second impurity dramatically changes the 
DOS obtained near a single impurity. In particular, we 
demonstrate that the number of resonance states in the 
DOS, as well as their frequency and spatial dependence 
change significantly as the distance between the impuri- 
ties or their orientation relative to the crystal lattice is 
varied. Since the latter effect arises from the momen- 



tum dependence of the superconducting gap, we argue 
that quantum interference is a novel tool to identify the 
symmetry of unconventional superconductors. This re- 
sult might be of particular importance for Sr2Ru04, an 
unconventional superconductor whose symmetry is still 
a topic of controversylld. While the study of quantum 
interference is not only of fundamental importance for 
our understanding of complex impurity structures, it can 
also clarify the origin of the resonances observed in the 
HTSC. In particular, we expect that the form of the res- 
onances arising from Kondo-screening of two magnetic 
impurities is different from those discussed below; work 
is currently under way to verify this conjecturellj. 

Starting—point for our calculations is the T-matrix 
formalisnxj which|-.we extended to treat scattering off 
multiple impuritiestrty in a d^2_y2-wa.ve superconductor. 



Quantum interference in-SiHiave superconductors was re- 
cently discussed in Refs.li3'L3. For simplicity we restrict 
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our considerations to two non-magnetic impurities; our 
formalism, however, allows the study of an arbitrary, but 
finite number of impurities. The study of more complex 
impurity structures, as well as that of piagnetic impu- 
rities, will be the focus of future workll3. Within the 
Nambu- formalism and for Matsubara frequencies, w„, 
the electronic Greens function in the presence of N im- 
purities is given by 

G{r,r',ujn) = G'o(r,r',cj„) 

N 

+ X! GQ{r,ri,uJn)T{ri,rj,uj„)Go{rj,r',uJn) , (1) 

where the T-matrix is obtained from the Bethe-Salpeter 
equation 



N 



+ Vri ^ Go (r,; , n , UJn)T(ri , Tj , UJn) 



(2) 



1=1 



For two non-magnetic impurities located at (i = 1,2) 
with Ar = |r2 — ri|, the scattering matrices are given by 
Vr^ = UiT^l2 with Ui being the non-magnetic scatter- 
ing potential and r the Pauli-matrices in Nambu-space. 
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The Greens function of the unperturbed (clean) system 
in momentum space is given by 



(3) 



For the electronic excitation spectrum in the normal state 
we take a form .that is characteristic of an optimally 
doped HTsdiaiilla and given by 

e\^ — —2t cos{kx) + cos{ky) — At' cos{kx) cos{ky) — ^ , 

(4) 

with t = 300 meV, t'/t — —0.4 as the nearest and next- 
nearest neighbor hopping integrals, respectively, and a 
chemical potential ^/t = —1.18, corresponding to 14% 
hole doping. Moreover, the superconducting gap with 
d2;2_j,2-symmetry iS|given by Ak = Aq [cos kx — cos ky] /2 
with Aq — 25 meVliZl. Our results presented below are 
qualitatively robust against changes in the form of et, 
or the size of Aq. We obtain the DOS, N{r,uj), from a 
numerical evaluation of Eqs.(0)- (HI) with A^(r,a;) = + 
A22 and Aii{r,uj) = -lmG^i{r,uJ + iS)/Tr. 

We briefly review some important features in the DOS 
near a single non-magnetic impurity in a (ij,2_y2-wave 
superconductorBi. The resulting diagonal T-matrix 
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11,22 — 



l-UoGo(r = 0,±u;) 



(5) 



where the upper (lower) sign applies to Tii(T22) and 
Go = [Go] 11 J exhibits a particle- {ujres < 0) and hole- 
like {ujres > 0) resonance. These resonances give rise 
to sharp peaks in the DOS only in the unitary limit 
{\ujres\/^o < 1) where U^^ = Re Go(0, icj^es)- 

In the presence of two impurities, the T-matrix, Eq.(^, 
is rather complex. However, in the limit _fo(Ar, lj) ^ 
Go(Ar, Lu), where Fq — [Go] 12 and for identical impurities 
with Ui^2 = Uq, the T-matrix simplifies considerably and 
is again diagonal. Defining 



S±iu) = {l-C/o[Go(0,c^)±Go(Ar,c^)]}- 



(6) 



one obtains (i j) 
2^11,22 (n,fj) 



Uo[S+i±Lo) + S-i±Lo)]/2 
f/oGo(Ar, ±Lu 



l-UoGo{0,±to) 



Tii,22{ri,n) (7) 



where the upper (lower) sign applies to Tii(T22). By 
comparing the T-matrices in Eqs.(^ and (0), we find that 
the presence of a second impurity splits the resonances 
of the single impurity case by J7oGo(Ar, w). Note that 
Go(Ar, w) does not only change with varying Ar, but 
also with the angle, a, between r2 — ri and the crystal 
X-axis, due to the momentum dependence of the super- 
conducting gap. Consequently, the energy and lifetime 
of the resonances depend on Ar and a. While all four 
5±-terms in Eq.(|7|) can possess resonances, those that are 
shifted to higher energies are highly overdamped and give 
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FIG. 1: DOS at R = (0, 0) as a function of Ar for two identi- 
cal impurities with Uo = 700 meV located at ri = (0, 0) and 
r2 — (Ar, 0). The lattice constant is set to ao = 1. 



rise only to oscillations in the DOS without the signature 
of a sharp peak. 

In what follows, we consider two identical impurities 
with scattering potential Ui^2 = Uq — 700 meV, coi; 
sponding to the unitary limit, in agreement with Refs.l 
While the specific form of the DOS changes with Uq, its 
qualitative features remain unchanged. To study the ef- 
fects of a and Ar on the DOS separately, we first consider 
for definiteness two impurities located along the crystal 
X-axis at ri = (0,0) and r2 = (Ar, 0) with a = 0. In 
Fig. ^ we present the DOS at R = (0,0), i.e., at one of 
the impurity sites, as a function of Ar. For comparison, 
we note that for a single impurity with Uq — 700 meV, 
the resonances are located ojres — ±1.5 meV. As Ar is 
varied, the DOS undergoes strong modifications. In par- 
ticular, the frequency of the resonances oscillates and at 
the same time, their energy width, or lifetime, changes. 
For a single impurity, the resonance frequency and width 
are correlated, such that as [t+iTinsl decreases, the fre- 
quency width decreases as welllira. In the case of two 
impurities, we find that \ujres \ and the lifetime of the res- 
onances are not necessarily correlated. For example, the 
resonance frequencies for both Ar = 2.0 and Ar = 3.5 
are LOres = ±4.0 meV, but the width of the resonances 
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FIG. 2: DOS as a function of spatial position R — {R, 0) for 
two impurities with Uo ~ 700 meV located at ri — (0, 0) and 
r2 = (Ar, 0) with (a) Ar = 2 and (b) Ar = 7. 



are considerably larger in the second case. Moreover, for 
some values of Ar, all resonances are very weak and, e.g., 
for Ar « 1 disappear almost completely. Note, that even 
for rather large values of Ar « 10, the DOS at R = (0, 0) 
is still affected by quantum interference and thus differ- 
ent from that obtained in the single impurity case. This 
result bears important implications for the interpretation 
of recent STM experimentsElEHO since it implies that the 
DOS near impurities in the two-dimensional HTSC can 
only be described within the single impurity framework 
if the impurity concentration is well below 1%. 

An additional important result of Fig. ^ is that the 
number of observable low-energy resonances changes with 
Ar. In particular, for Ar < 6, only two sharp low-energy 
resonances can be clearly identified. This effect becomes 
particularly evident when one considers the spatial de- 
pendence of the DOS for fixed Ar, as shown in Fig. ||. 
Here, we plot the DOS as a function of R = {R, 0) for two 
impurities located at ri — (0,0) and r2 — (Ar, 0). The 
uppermost curve corresponds to the midpoint between 
the two impurities, the dashed line represents the DOS at 
r2 . For Ar = 2 (Fig. ||a) , there exist only two low-energy 
resonances at ujres i 4 meV. In contrast, for Ar = 7 
(Fig. ||b), we find two broader resonances at ujres — ±3 
meV, and two sharper resonances at cOres = ±0.25 meV. 
Note, that the resonances for Ar — 2 a,t ujres^4: meV have 
a considerably larger amplitude than those for Ar = 7 at 
LOres = i3 meV. This result again differs from the sin- 
gle impurity case, where the resonance with the smaller 
{ujresl always possesses a larger amplitude in the DOS. 

While sharp resonances can only be identified for 
ujres <C Ao, oscillations in the DOS exist for basically 
all energies \uj\ < Ag. To study these oscillations in more 
detail, we present in Fig. || the DOS along R = {R, 0) for 
various frequencies and the same impurity arrangement 
as in Fig. ^a; the locations of the impurities are indicated 
by arrows. The solid and dashed lines in Fig. Ha represent 
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FIG. 3: (DOS along R = (R, 0) for the same impurity ar- 
rangement as in Fig. ^a. The positions of the impurities 
are indicated by arrows, (a) DOS at u>res ~ ±4 meV. The 
open and filled squares present the DOS at tJres = ±4 meV 
along the lattice diagonal with R — {R, R). Inset: DOS along 
R = {R, 0) for a single impurity at its resonance frequency 
LUres = ±1.5 meV. (b) DOS at frequencies \u!res\ < |w| < Aq. 



the DOS at ojres — ±4 meV, corresponding to peak (1) 
and (2) in Fig. ||a. The DOS at \ujres\ but along the lat- 
tice diagonal, i.e., for R — (R, R), is shown as open and 
filled squares. Similar to the single impurity case, the 
amplitude of the resonances is much weaker along the 
direction of the superconducting gap nodes, than along 
the anti-nodal direction. The inset shows the DOS along 
R = (i?,0) for a single impurity with Uq = 700 meV 
and resonance energy Ures = ±1.5 meV. A comparison 
shows that the amplitude of the DOS oscillations induced 
by two impurities is much larger than those induced by 
a single impurity (the overall scale in the inset is three 
times smaller than in the main figure). Moreover, in the 
two impurity case, the DOS exhibits significant oscilla- 
tions at much larger distances from the impurities than 
in the single impurity case. This is particularly evident 
when comparing the particle-like resonances, where in 
the two impurity case, the amplitude of the oscillations 
is still large at a distance to the nearest impurity, r„, of 
about r„ « 4—5, while in the single-impurity case, the os- 
cillations are already substantially reduced at r„ w 2. In 
Fig. 1^ we present the DOS along R — {R, 0) for several 
frequencies with Iwresl < |<^| < Ag. While there exist no 
evidence for a resonance at higher energies, we still find 
considerable oscillations in the DOS. As |a;| decreases, 
the wave-vectors of these oscillations decreases, as can 
clearly be seen from the shift of the peaks around R = 3 
and 4. Thus, the DOS oscillations exhibit a dispersion, 
similar to the results obtained in Rcf.Ll. 

Due to the momentum dependence of the supercon- 
ducting d^2_j,2-gap, the DOS changes when the orien- 
tation of the two impurities with respect to the crystal 
lattice is varied. In particular, since the gap vanishes 
along the lattice diagonal, we expect the largest devia- 



4 




-30 -25 -20 -15 -10 -5 5 10 15 20 25 30 

Frequency [meV] 

FIG. 4: DOS at (a) R = (0,0) and (b) R = (1,0) for two 
impurities with U — 700 meV. One impurity is located at 
ri = (0, 0), the other one either at r2 — (3, 3) (dotted line) or 
at r2 = (4,0) (dashed line). 



tions from the results for a = shown in Figs. |l|-^ when 
the impurities are located along the lattice diagonal, i.e., 
for a — tt/4:. To study the changes in the DOS due to 
variations in a and to eliminate effects due to a varying 



Ar, we chose two different impurity arrangements which 
can be realized experimentally, and possess almost iden- 
tical values for Ar. In the first case, the impurities are 
located along the lattice diagonal (a = 7r/4) at ri = (0, 0) 
and r2 = (3,3) (dotted line, Ar « 4.2). In the second 
arrangement, the impurities are aligned along the crystal 
X-axis {a = 0) at ri = (0, 0) and r2 = (4, 0) (dashed line, 
Ar = 4). We verified that the DOS for rs = (4.2,0), 
which would yield identical values of Ar, but is experi- 
mentally not realizable, differs only slightly from that for 
r2 = (4,0). The DOS at R = (0,0) and R = (1,0) is 
shown in Fig. ^. While the DOS for a — tt/A possesses 
three distinct resonances, only two resonance peaks are 
observable for a = 0. Moreover, for a = Q, the reso- 
nance states are located at higher frequencies and are 
much broader. The qualitative differences in the DOS 
between these two different impurity arrangements can 
be directly traced back to the vanishing of F(Ar, oj) for 
a = 7r/4, and its finite, complex value for a — 0. Thus, 
the symmetry of the superconducting gap is directly re- 
flected in the changes which the DOS undergoes when the 
orientation of the impurities relative to the crystal lattice 
is varied. This dependence provides a new tool to identify 
the symmetry of unconventional superconductors. 

In summary, we studied quantum interference of elec- 
tronic waves that are scattered by two non-magnetic im- 
purities in a d3,2_y2-superconductor. We show that the 
number of resonance states in the DOS, as well as their 
frequency and spatial dependence changes significantly as 
the distance between the impurities or their orientation 
relative to the crystal lattice is varied. The latter result 
provides a novel tool to identify the symmetry of uncon- 
ventional superconductors, such as Sr2Ru04, where the 
symmetry of the superconducting state is still a topic of 
controversy. 

We would like to thank J.C. Davis and A. de Lozanne 
for stimulating discussions. 
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